RAB32 Ser71Arg in autosomal dominant Parkinson's disease: linkage, association, and functional analyses.

BACKGROUND: Parkinson's disease is a progressive neurodegenerative disorder with multifactorial causes, among which genetic risk factors play a part. The RAB GTPases are regulators and substrates of LRRK2, and variants in the LRRK2 gene are important risk factors for Parkinson's disease. We aimed to explore genetic variability in RAB GTPases within cases of familial Parkinson's disease. METHODS: We did whole-exome sequencing in probands from families in Canada and Tunisia with Parkinson's disease without a genetic cause, who were recruited from the Centre for Applied Neurogenetics (Vancouver, BC, Canada), an international consortium that includes people with Parkinson's disease from 36 sites in 24 countries. 61 RAB GTPases were genetically screened, and candidate variants were genotyped in relatives of the probands to assess disease segregation by linkage analysis. Genotyping was also done to assess variant frequencies in individuals with idiopathic Parkinson's disease and controls, matched for age and sex, who were also from the Centre for Applied Neurogenetics but unrelated to the probands or each other. All participants were aged 18 years or older. The sequencing and genotyping findings were validated by case-control association analyses using bioinformatic data obtained from publicly available clinicogenomic databases (AMP-PD, GP2, and 100 000 Genomes Project) and a private German clinical diagnostic database (University of Tübingen). Clinical and pathological findings were summarised and haplotypes were determined. In-vitro studies were done to investigate protein interactions and enzyme activities. FINDINGS: Between June 1, 2010, and May 31, 2017, 130 probands from Canada and Tunisia (47 [36%] female and 83 [64%] male; mean age 72·7 years [SD 11·7; range 38-96]; 109 White European ancestry, 18 north African, two east Asian, and one Hispanic] underwent whole-exome sequencing. 15 variants in RAB GTPase genes were identified, of which the RAB32 variant c.213C>G (Ser71Arg) cosegregated with autosomal dominant Parkinson's disease in three families (nine affected individuals; non-parametric linkage Z score=1·95; p=0·03). 2604 unrelated individuals with Parkinson's disease and 344 matched controls were additionally genotyped, and five more people originating from five countries (Canada, Italy, Poland, Turkey, and Tunisia) were identified with the RAB32 variant. From the database searches, in which 6043 individuals with Parkinson's disease and 62 549 controls were included, another eight individuals were identified with the RAB32 variant from four countries (Canada, Germany, UK, and USA). Overall, the association of RAB32 c.213C>G (Ser71Arg) with Parkinson's disease was significant (odds ratio [OR] 13·17, 95% CI 2·15-87·23; p=0·0055; I2=99·96%). In the people who had the variant, Parkinson's disease presented at age 54·6 years (SD 12·75, range 31-81, n=16), and two-thirds had a family history of parkinsonism. RAB32 Ser71Arg heterozygotes shared a common haplotype, although penetrance was incomplete. Findings in one individual at autopsy showed sparse neurofibrillary tangle pathology in the midbrain and thalamus, without Lewy body pathology. In functional studies, RAB32 Arg71 activated LRRK2 kinase to a level greater than RAB32 Ser71. INTERPRETATION: RAB32 Ser71Arg is a novel genetic risk factor for Parkinson's disease, with reduced penetrance. The variant was found in individuals with Parkinson's disease from multiple ethnic groups, with the same haplotype. In-vitro assays show that RAB32 Arg71 activates LRRK2 kinase, which indicates that genetically distinct causes of familial parkinsonism share the same mechanism. The discovery of RAB32 Ser71Arg also suggests several genetically inherited causes of Parkinson's disease originated to control intracellular immunity. This shared aetiology should be considered in future translational research, while the global epidemiology of RAB32 Ser71Arg needs to be assessed to inform genetic counselling. FUNDING: National Institutes of Health, the Canada Excellence Research Chairs program, Aligning Science Across Parkinson's, the Michael J Fox Foundation for Parkinson's Research, and the UK Medical Research Council.

A pathogenic variant in RAB32 causes autosomal dominant Parkinson's disease and activates LRRK2 kinase.

Background: Parkinson's disease (PD) is a progressive neurodegenerative disorder. Mendelian forms have revealed multiple genes, with a notable emphasis on membrane trafficking; RAB GTPases play an important role in PD as a subset are both regulators and substrates of LRRK2 protein kinase. To explore the role of RAB GTPases in PD, we undertook a comprehensive examination of their genetic variability in familial PD. Methods: Affected probands from 130 multi-incident PD families underwent whole-exome sequencing and genotyping, Potential pathogenic variants in 61 RAB GTPases were genotyped in relatives to assess disease segregation. These variants were also genotyped in a larger case-control series, totaling 3,078 individuals (2,734 with PD). The single most significant finding was subsequently validated within genetic data (6,043 with PD). Clinical and pathologic findings were summarized for gene-identified patients, and haplotypes were constructed. In parallel, wild-type and mutant RAB GTPase structural variation, protein interactions, and resultant enzyme activities were assessed. Findings: We found RAB32 c.213C>G (Ser71Arg) to co-segregate with autosomal dominant parkinsonism in three multi-incident families. RAB32 Ser71Arg was also significantly associated with PD in case-control samples: genotyping and database searches identified thirteen more patients with the same variant that was absent in unaffected controls. Notably, RAB32 Ser71Arg heterozygotes share a common haplotype. At autopsy, one patient had sparse neurofibrillary tangle pathology in the midbrain and thalamus, without Lewy body pathology. In transfected cells the RAB32 Arg71 was twice as potent as Ser71 wild type to activate LRRK2 kinase. Interpretation: Our study provides unequivocal evidence to implicate RAB32 Ser71Arg in PD. Functional analysis demonstrates LRRK2 kinase activation. We provide a mechanistic explanation to expand and unify the etiopathogenesis of monogenic PD. Funding: National Institutes of Health, the Canada Excellence Research Chairs program, Aligning Science Across Parkinson's, the Michael J. Fox Foundation for Parkinson's Research, and the UK Medical Research Council.

Mitochondrial DNA heteroplasmy distinguishes disease manifestation in PINK1/PRKN-linked Parkinson's disease.

Biallelic mutations in PINK1/PRKN cause recessive Parkinson's disease. Given the established role of PINK1/Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA integrity and inflammation as disease modifiers in carriers of mutations in these genes. Mitochondrial DNA integrity was investigated in a large collection of biallelic (n = 84) and monoallelic (n = 170) carriers of PINK1/PRKN mutations, idiopathic Parkinson's disease patients (n = 67) and controls (n = 90). In addition, we studied global gene expression and serum cytokine levels in a subset. Affected and unaffected PINK1/PRKN monoallelic mutation carriers can be distinguished by heteroplasmic mitochondrial DNA variant load (area under the curve = 0.83, CI 0.74-0.93). Biallelic PINK1/PRKN mutation carriers harbour more heteroplasmic mitochondrial DNA variants in blood (P = 0.0006, Z = 3.63) compared to monoallelic mutation carriers. This enrichment was confirmed in induced pluripotent stem cell-derived (controls, n = 3; biallelic PRKN mutation carriers, n = 4) and post-mortem (control, n = 1; biallelic PRKN mutation carrier, n = 1) midbrain neurons. Last, the heteroplasmic mitochondrial DNA variant load correlated with IL6 levels in PINK1/PRKN mutation carriers (r = 0.57, P = 0.0074). PINK1/PRKN mutations predispose individuals to mitochondrial DNA variant accumulation in a dose- and disease-dependent manner.

A Whole-Genome Sequencing Study Implicates GRAMD1B in Multiple Sclerosis Susceptibility.

While the role of common genetic variants in multiple sclerosis (MS) has been elucidated in large genome-wide association studies, the contribution of rare variants to the disease remains unclear. Herein, a whole-genome sequencing study in four affected and four healthy relatives of a consanguineous Italian family identified a novel missense c.1801T > C (p.S601P) variant in the GRAMD1B gene that is shared within MS cases and resides under a linkage peak (LOD: 2.194). Sequencing GRAMD1B in 91 familial MS cases revealed two additional rare missense and two splice-site variants, two of which (rs755488531 and rs769527838) were not found in 1000 Italian healthy controls. Functional studies demonstrated that GRAMD1B, a gene with unknown function in the central nervous system (CNS), is expressed by several cell types, including astrocytes, microglia and neurons as well as by peripheral monocytes and macrophages. Notably, GRAMD1B was downregulated in vessel-associated astrocytes of active MS lesions in autopsied brains and by inflammatory stimuli in peripheral monocytes, suggesting a possible role in the modulation of inflammatory response and disease pathophysiology.

Association of Essential Tremor With Novel Risk Loci: A Genome-Wide Association Study and Meta-analysis.

Importance: Essential tremor (ET) is one of the most common movement disorders, affecting 5% of the general population older than 65 years. Common variants are thought to contribute toward susceptibility to ET, but no variants have been robustly identified. Objective: To identify common genetic factors associated with risk of ET. Design, Setting, and Participants: Case-control genome-wide association study. Inverse-variance meta-analysis was used to combine cohorts. Multicenter samples collected from European populations were collected from January 2010 to September 2019 as part of an ongoing study. Included patients were clinically diagnosed with or reported having ET. Control individuals were not diagnosed with or reported to have ET. Of 485 250 individuals, data for 483 054 passed data quality control and were used. Main Outcomes and Measures: Genotypes of common variants associated with risk of ET. Results: Of the 483 054 individuals included, there were 7177 with ET (3693 [51.46%] female; mean [SD] age, 62.66 [15.12] years), and 475 877 control individuals (253 785 [53.33%] female; mean [SD] age, 56.40 [17.6] years). Five independent genome-wide significant loci and were identified and were associated with approximately 18% of ET heritability. Functional analyses found significant enrichment in the cerebellar hemisphere, cerebellum, and axonogenesis pathways. Genetic correlation (r), which measures the degree of genetic overlap, revealed significant common variant overlap with Parkinson disease (r, 0.28; P = 2.38 × 10-8) and depression (r, 0.12; P = 9.78 × 10-4). A separate fine-mapping of transcriptome-wide association hits identified genes such as BACE2, LRRN2, DHRS13, and LINC00323 in disease-relevant brain regions, such as the cerebellum. Conclusions and Relevance: The results of this genome-wide association study suggest that a portion of ET heritability can be explained by common genetic variation and can help identify new common genetic risk factors for ET.

Elevated levels of serum CD5 antigen-like protein distinguish secondary progressive multiple sclerosis from other disease subtypes.

CD5 antigen-like (CD5L) protein is a macrophage-secreted protein with roles in immunomodulation and lipid homeostasis. We compared serum CD5L levels in healthy controls to individuals diagnosed with clinically isolated syndrome, relapsing remitting (RR), secondary progressive (SP), and primary progressive (PP) multiple sclerosis (MS). CD5L was increased in SPMS relative to controls, RRMS, and PPMS. SPMS CD5L was associated with longer disease duration independent of age, sex, or disease severity. The positive relationship between CD5L and disease duration in SPMS suggests a chronic peripheral inflammatory profile compared to other subtypes, particularly PPMS, warranting investigation of functional roles for CD5L in MS.

Regulator of calcineurin 1 is a novel RNA-binding protein to regulate neuronal apoptosis.

Posttranscriptional regulation of gene expression plays an important role in the maturation, transport, stability and translation of coding and noncoding RNAs. RNA-binding protein (RBP) is a key factor of the regulation. Regulator of calcineurin 1 (RCAN1) is a multifunctional protein involved in neurodegeneration, mitochondrial dysfunction, inflammation and protein glycosylation, and plays an important role in the pathogenesis of Down syndrome and Alzheimer's disease. In this report, we discovered that RCNA1 is a novel RNA-binding protein. A 23 nucleotide sequence of adenine nucleotide translocator (ANT1) mRNA was identified as the binding motif of RCAN1. Furthermore, we found that R1SR13, as the RNA aptamer of RCAN1 identified by SELEX, blocked RCAN1-induced inhibition of the nuclear factor of activated T cells (NFAT) and NF-κB signaling pathways, and reduced neuronal apoptosis. Taken together, our results demonstrate that RCAN1 is a novel RNA-binding protein and the RNA aptamer of RCAN1 plays a neuroprotective role.

Exome-wide rare variant analysis in familial essential tremor.

INTRODUCTION: Essential tremor (ET) is one of the most common movement disorders. Despite its high prevalence and heritability, its genetic etiology remains elusive with only a few susceptibility genes identified and poorly replicated. Our aim was to find novel candidate genes involved in ET predisposition through whole exome sequencing. METHODS: We studied eight multigenerational families (N = 40 individuals) with an autosomal-dominant inheritance using a comprehensive strategy combining whole exome sequencing followed by case-control association testing of prioritized variants in a separate cohort comprising 521 ET cases and 596 controls. We further performed gene-based burden analyses in an additional dataset comprising 789 ET patients and 770 healthy individuals to investigate whether there was an enrichment of rare deleterious variants within our candidate genes. RESULTS: Fifteen variants co-segregated with disease status in at least one of the families, among which rs749875462 in CCDC183, rs535864157 in MMP10 and rs114285050 in GPR151 showed a nominal association with ET. However, we found no significant enrichment of rare variants within these genes in cases compared with controls. Interestingly, MMP10 protein is involved in the inflammatory response to neuronal damage and has been previously associated with other neurological disorders. CONCLUSIONS: We prioritized a set of promising genes, especially MMP10, for further genetic and functional studies in ET. Our study suggests that rare deleterious coding variants that markedly increase susceptibility to ET are likely to be found in many genes. Future studies are needed to replicate and further infer biological mechanisms and potential disease causality for our identified genes.

Enhanced expression of complement and microglial-specific genes prior to clinical progression in the MOG-experimental autoimmune encephalomyelitis model of multiple sclerosis.

Understanding the biological changes responsible for failures in repair and the development of progressive MS is paramount for therapeutic intervention. In a well characterized experimental autoimmune encephalomyelitis (EAE) model of MS the clinical phenotype features an acute attack with partial recovery followed by a chronic or progressive disease phase. Neuropathology-focused gene expression profiles were generated from spinal cord, hindbrain and forebrain of mice 25 days after the induction of EAE, the time when recovery plateaus and transitions to a chronic or worsening phase. Differences in gene expression were most pronounced in the spinal cord of EAE mice compared to sham-immunized animals, with a subset of genes also found to be differentially expressed in the hindbrain and the forebrain, albeit with smaller fold-changes in expression. Our data suggests that changes in complement components, chemoattractant cytokines and especially enrichment in microglial cells may be the primary drivers of processes that limit recovery in EAE.

Genetic analysis of nucleotide-binding leucine-rich repeat (NLR) receptors in multiple sclerosis.

Genetic and functional analyses of the inflammasome suggest a role for this multiprotein complex in the biological mechanisms leading to the onset and progression of multiple sclerosis (MS). Nucleotide-binding, leucine-rich repeat (NLR) receptors trigger the activation and assembly of specific inflammasomes in response to danger signals. Mining exome sequencing data from 326 MS patients identified 17 rare missense or nonsense variants in NLR family pyrin domain containing 1 (NLRP1), NLRP3, NLRP6, NLRP7 and NLR family CARD domain containing 4 (NLRC4). Genotyping these variants in 2503 MS cases and 1076 healthy controls did not result in statistically significant differences between groups, and segregation analysis within MS families was largely unsupportive of co-segregation of these variants with disease. However, the identification of MS patients harboring rare homozygote variants in NLRP1 (p.Ile601Phe and p.Ser1387Ile), a variant in NLRP3 (p.Leu832Ile) resulting in the substitution of a critical amino acid for the formation of its leucine-rich repeat domain, and several MS patients with NLRC4 variants (p.Arg310Ter and p.Glu600Ter) causing protein truncations suggest that rare protein-altering variants in inflammasome-activating NLR receptors may contribute to MS risk.

Oligodendrocyte ARNT2 expression is altered in models of MS.

OBJECTIVE: We examined expression of aryl hydrocarbon receptor nuclear translocator 2 (ARNT2), a basic-loop-helix transcription factor implicated in neuronal development and axonal health, in oligodendrocyte (OL) cultures and over the course of chronic experimental autoimmune encephalomyelitis (EAE), the murine model of multiple sclerosis (MS). METHODS: We assessed OL ARNT2 expression in EAE compared with sham-immunized controls and also in OL primary cultures and over the course of dibutyryl cyclic adenosine monophosphate (dbcAMP)-mediated maturation of the immortalized Oli-neu cell line. We also tested the functional role of ARNT2 in influencing OL characteristics using small interfering RNA (siRNA). RESULTS: ARNT2 is localized to Olig2+ cells in healthy spinal cord gray and white matter. Despite a significant expansion of Olig2+ cells in the white matter at peak disease, ARNT2 is reduced by almost half in OLs, along with a reduction in the percentage of ARNT2+/Olig2+ cells. Mature OLs in mixed cortical cultures or OLs matured from embryonic progenitors express negligible ARNT2. Similarly, Oli-neu cells express high levels of ARNT2, which are reduced following dbcAMP maturation. siRNA-mediated knockdown of ARNT2 affected OL viability, which led to an enrichment of myelin-producing OLs. CONCLUSION: The analysis of ARNT2 expression in OLs demonstrates that OL ARNT2 expression is altered in EAE and during OL maturation. Findings point to ARNT2 as an important mediator of OL viability and differentiation and warrant further characterization as a target for intervention in demyelinating disorders such as MS.

NLRX1 inhibits the early stages of CNS inflammation and prevents the onset of spontaneous autoimmunity.

Nucleotide-binding, leucine-rich repeat containing X1 (NLRX1) is a mitochondria-located innate immune sensor that inhibits major pro-inflammatory pathways such as type I interferon and nuclear factor-κB signaling. We generated a novel, spontaneous, and rapidly progressing mouse model of multiple sclerosis (MS) by crossing myelin-specific T-cell receptor (TCR) transgenic mice with Nlrx1-/- mice. About half of the resulting progeny developed spontaneous experimental autoimmune encephalomyelitis (spEAE), which was associated with severe demyelination and inflammation in the central nervous system (CNS). Using lymphocyte-deficient mice and a series of adoptive transfer experiments, we demonstrate that genetic susceptibility to EAE lies within the innate immune compartment. We show that NLRX1 inhibits the subclinical stages of microglial activation and prevents the generation of neurotoxic astrocytes that induce neuronal and oligodendrocyte death in vitro. Moreover, we discovered several mutations within NLRX1 that run in MS-affected families. In summary, our findings highlight the importance of NLRX1 in controlling the early stages of CNS inflammation and preventing the onset of spontaneous autoimmunity.

Exome sequencing in multiple sclerosis families identifies 12 candidate genes and nominates biological pathways for the genesis of disease.

Multiple sclerosis (MS) is an inflammatory disease of the central nervous system characterized by myelin loss and neuronal dysfunction. Although the majority of patients do not present familial aggregation, Mendelian forms have been described. We performed whole-exome sequencing analysis in 132 patients from 34 multi-incident families, which nominated likely pathogenic variants for MS in 12 genes of the innate immune system that regulate the transcription and activation of inflammatory mediators. Rare missense or nonsense variants were identified in genes of the fibrinolysis and complement pathways (PLAU, MASP1, C2), inflammasome assembly (NLRP12), Wnt signaling (UBR2, CTNNA3, NFATC2, RNF213), nuclear receptor complexes (NCOA3), and cation channels and exchangers (KCNG4, SLC24A6, SLC8B1). These genes suggest a disruption of interconnected immunological and pro-inflammatory pathways as the initial event in the pathophysiology of familial MS, and provide the molecular and biological rationale for the chronic inflammation, demyelination and neurodegeneration observed in MS patients.

Absence of Mutation Enrichment for Genes Phylogenetically Conserved in the Olivocerebellar Motor Circuitry in a Cohort of Canadian Essential Tremor Cases.

Essential Tremor is a prevalent neurological disorder of unknown etiology. Studies suggest that genetic factors contribute to this pathology. To date, no causative mutations in a gene have been reproducibly reported. All three structures of the olivocerebellar motor circuitry have been linked to Essential Tremor. We postulated that genes enriched for their expression in the olivocerebellar circuitry would be more susceptible to harbor mutations in Essential Tremor patients. A list of 11 candidate genes, enriched for their expression in the olivocerebellar circuitry, was assessed for their variation spectrum and frequency in a cohort of Canadian Essential Tremor cases. Our results from this list of 11 candidate genes do not support an association for Essential Tremor in our cohort of Canadian cases. The heterogenic nature of ET and modest size of the cohort used in this study are two confounding factors that could explain these results.

The Role of Rare Coding Variants in Parkinson's Disease GWAS Loci.

Introduction: Genome-wide association studies (GWAS) have identified multiple loci associated with Parkinson's disease (PD) risk. The presence of rare variants within these loci that may account for the increased susceptibility requires further investigation. Methods: Using exome sequencing, we performed a comprehensive rare variant screen of genes located within 56 novel PD loci. We first analyzed exomes from 109 subjects in the discovery cohort (85 diagnosed with PD and 24 healthy controls) and filtered for rare coding variants with minor allele frequency <0.01 and combined annotation-dependent depletion > 15. Further investigation of exome data from a replication cohort of 2,859 European patients with PD (International Parkinson's Disease Genomics Consortium) and 24,146 non-Finnish European controls from gnomAD were used for association testing of specific rare variants found in the discovery cohort. Results: Our genetic screening identified 54 potential disease-relevant variants in 71 genes in 109 subjects. Six out of 54 variants were found in two or more patients and were not observed in healthy controls: DNAH1 p.A3639T, STAB1 p.S1089G, ANK2 p.V3634D, ANK2 p.R3906W, SH3GL2 p.G276V, and NOD2 p.G908R. Replication in the International Parkinson's Disease Genomics Consortium (IPDGC) confirmed the association with PD risk for three out of the six identified variants (STAB1 p.S1089G, SH3GL2 p.G276V, and NOD2 p.G908R) (p < 10-3). Conclusion: Our study suggests that some of the associations identified in PD risk loci can be ascribed to rare variants with likely functional effects that modify PD risk.

Whole-Exome Sequencing of an Exceptional Longevity Cohort.

Centenarians represent a unique cohort to study the genetic basis for longevity and factors determining the risk of neurodegenerative disorders, including Alzheimer's disease (AD). The estimated genetic contribution to longevity is highest in centenarians and super-cententenarians, but few genetic variants have been shown to clearly impact this phenotype. While the genetic risk for AD and other dementias is now well understood, the frequency of known dementia risk variants in centenarians is not fully characterized. To address these questions, we performed whole-exome sequencing on 100 individuals of 98-108 years age in search of genes with large effect sizes towards the exceptional aging phenotype. Overall, we were unable to identify a rare protein-altering variant or individual genes with an increased burden of rare variants associated with exceptional longevity. Gene burden analysis revealed three genes of nominal statistical significance associated with extreme aging, including LYST, MDN1, and RBMXL1. Several genes with variants conferring an increased risk for AD and other dementias were identified, including TREM2, EPHA1, ABCA7, PLD3, MAPT, and NOTCH3. Larger centenarian studies will be required to further elucidate the genetic basis for longevity, and factors conferring protection against age-dependent neurodegenerative syndromes.

TPP2 mutation associated with sterile brain inflammation mimicking MS.

OBJECTIVE: To ascertain the genetic cause of a consanguineous family from Syria suffering from a sterile brain inflammation mimicking a mild nonprogressive form of MS. METHODS: We used homozygosity mapping and next-generation sequencing to detect the disease-causing gene in the affected siblings. In addition, we performed RNA and protein expression studies, enzymatic activity assays, immunohistochemistry, and targeted sequencing of further MS cases from Austria, Germany, Canada and Jordan. RESULTS: In this study, we describe the identification of a homozygous missense mutation (c.82T>G, p.Cys28Gly) in the tripeptidyl peptidase II (TPP2) gene in all 3 affected siblings of the family. Sequencing of all TPP2-coding exons in 826 MS cases identified one further homozygous missense variant (c.2027C>T, p.Thr676Ile) in a Jordanian MS patient. TPP2 protein expression in whole blood was reduced in the affected siblings. In contrast, TPP2 protein expression in postmortem brain tissue from MS patients without TPP2 mutations was highly upregulated. CONCLUSIONS: The homozygous TPP2 mutation (p.Cys28Gly) is likely responsible for the inflammation phenotype in this family. TPP2 is an ubiquitously expressed serine peptidase that removes tripeptides from the N-terminal end of longer peptides. TPP2 is involved in various biological processes including the destruction of major histocompatibility complex Class I epitopes. Recessive loss-of-function mutations in TPP2 were described in patients with Evans syndrome, a rare autoimmune disease affecting the hematopoietic system. Based on the gene expression results in our MS autopsy brain samples, we further suggest that TPP2 may play a broader role in the inflammatory process in MS.

Exome sequencing study in patients with multiple sclerosis reveals variants associated with disease course.

BACKGROUND: It remains unclear whether disease course in multiple sclerosis (MS) is influenced by genetic polymorphisms. Here, we aimed to identify genetic variants associated with benign and aggressive disease courses in MS patients. METHODS: MS patients were classified into benign and aggressive phenotypes according to clinical criteria. We performed exome sequencing in a discovery cohort, which included 20 MS patients, 10 with benign and 10 with aggressive disease course, and genotyping in 2 independent validation cohorts. The first validation cohort encompassed 194 MS patients, 107 with benign and 87 with aggressive phenotypes. The second validation cohort comprised 257 patients, of whom 224 patients had benign phenotypes and 33 aggressive disease courses. Brain immunohistochemistries were performed using disease course associated genes antibodies. RESULTS: By means of single-nucleotide polymorphism (SNP) detection and comparison of allele frequencies between patients with benign and aggressive phenotypes, a total of 16 SNPs were selected for validation from the exome sequencing data in the discovery cohort. Meta-analysis of genotyping results in two validation cohorts revealed two polymorphisms, rs28469012 and rs10894768, significantly associated with disease course. SNP rs28469012 is located in CPXM2 (carboxypeptidase X, M14 family, member 2) and was associated with aggressive disease course (uncorrected p value < 0.05). SNP rs10894768, which is positioned in IGSF9B (immunoglobulin superfamily member 9B) was associated with benign phenotype (uncorrected p value < 0.05). In addition, a trend for association with benign phenotype was observed for a third SNP, rs10423927, in NLRP9 (NLR family pyrin domain containing 9). Brain immunohistochemistries in chronic active lesions from MS patients revealed expression of IGSF9B in astrocytes and macrophages/microglial cells, and expression of CPXM2 and NLRP9 restricted to brain macrophages/microglia. CONCLUSIONS: Genetic variants located in CPXM2, IGSF9B, and NLRP9 have the potential to modulate disease course in MS patients and may be used as disease activity biomarkers to identify patients with divergent disease courses. Altogether, the reported results from this study support the influence of genetic factors in MS disease course and may help to better understand the complex molecular mechanisms underlying disease pathogenesis.